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Abstract 

In this article, we prove decorrelation estimates for the eigenvalues of 
a ID discrete tight binding model near two distinct energies in the 
localized regime. Consequently, with an arbitrary, fixed number n, 
the asymptotic independence for local level statistics near n distinct 
energies is obtained. 



1 Introduction 

The present paper deals with the following lattice Hamiltonian with off- 
diagonal disorder in dimension 1: for u = {u{n)}n<=z ^ ^'^i'^), set 

{H^u){n) = iOniu{n) - u{n + 1)) - ujn~i{u{n - 1) - u{n)). (1.1) 

This model appears in the description of waves (light, acoustic waves, etc) 
which propagate through a disordered, discrete medium (c.f. [1]). 
We can see {cj„}„gz in this model as weights of bonds of the lattice Z. 
Throughout this article, we assume that u := {Un}nez are non- negative i.i.d. 
random variables (r.v.'s for short) with a bounded, compactly supported 
density p. 

In addition, from Section 1 to Section 3, we assume more that Un G [ao; f^o] 
for all n G Z where /3o > ao > 0. In Section 4, we will comment on relaxing 
the hypothesis of the lower bound of r.v.'s u. 
It is known that (see [ti]): 

• the operator H^^ admits an almost sure spectrum S := a^H^^) = [0, 4/3o]. 

• has an integrated density of states defined as follows: 
cj— a.s., the following limit exists and is u independent: 

, {e.v. of iJ^ less than E} ^ „ . , 

N{E) := lim ^ ^ for a.e. E. (1.2) 

|A|->-+oo A 



1 



1 INTRODUCTION 



2 



As a direct consequence of the Wegner estimate (see Theorem 2.1 in 
Section 2), N{E) is defined everywhere in M and absolutely continuous 
w.r.t. Lebesgue measure with a bounded derivative z^(-E') called the 
density of states of H^. 

In the present paper, we follow a usual way to study various statistics related 
to random operators. We restrict the operator iJ^ on some interval A C Z 
of finite length with some boundary condition and obtain a finite-volume 
operator which is denoted by if^(A). Then, we study diverse statistics for 
this operator in the limit when |A| goes to infinity. 

Throughout this paper, the boundary condition to define H^{A) is always 
the periodic boundary condition. For example, if A = [1,A^], the operator 
Hi^{A) is a self-adjoint N x N matrix of the following form: 

/ un + ooi —uJi . . . —ujn \ 

—UJ\ + ijJ2 —002 ■ ■ ■ 

... IjJn-2 + ojn-i —ojn-1 

y -UJN ... -W7v_i UJn^I+OJn) 

For L G N, let A = A^ := L] be a large interval in Z and |A| := (2L + 1) 
be its cardinality. 

We will denote the eigenvalues of H^{K) ordered increasingly and repeated 
according to multiplicity by i?i(a;. A) ^ £'2(w. A) ^ ■ ■ ■ ^ E\^{uj,K). 
Let / be the localized region in S where the finite-volume fractional- moment 
criteria for localization are satisfied for the finite- volume operators Hi_j{A) 
when |A| large enough (see Section 2 for more details). In this region, the 
spectrum of H^^ is pure point and the corresponding eigenfunctions are ex- 
ponentially decay. 

Pick E a positive energy in I with ^{E) > 0, and define the local level 
statistics near E as follows 

|A| 

E{^,E,u,A) = J2^UE.^,m) (1-3) 

n=l 

where 

UE, CO, A) = \A\u{E){Er,iuj, A)-E). (1.4) 

For the model (1.1), It is known that the weak limit of the above point 
process is a Poisson point process: 

Theorem 1.1. [[6], [8]] Assume that E is a positive energy in I with ^{E) > 
0. 

Then, when \ A\ — )■ +oo, the point process E{^,E,uj,A) converges weakly to a 
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Poisson point process with the intensity 1 i.e., for {Uj)i<^j<^j, Uj C M bounded 
measurable and Uji fl f/j = if j ^ j' and {kj)i<^j<^j G N"^, we have 



lim 

|Al^+oo 



P 



\(i^{j;^j{E,io,A)eU,} =kjj j 



krV^ 






-n 


kj\] 


i=i 



k,\ 



Recently, for the ID discrete Anderson model, Klopp [9] showed more that 
if we pick two fixed, distinct energies E and E' in the localized regime, their 
two corresponding point processes 'E[^,E,uj, K) and w. A) converge 

weakly, respectively to two independent Poisson point processes. In other 
words, the limits of w. A) and ii^', A) are stochastically inde- 

pendent. 

From the argument in [9], the above statement holds true if one can prove a 
so-called decorrelation estimate. 

That is exactly what we want to carry out here for the ID discrete lattice 
Hamiltonian with off-diagonal disorder (1.1). Our decorrelation estimate is 
the following: 

Theorem 1.2. Let E,E' be two positive, distinct energies in the localized 
regime. Pick (5 G (1/2, 1) and a G (0, 1) . Then, for any c > 0, there exists 
C > such that, for L large enough and cV^ ^ / ^ L'^ /c, one has 

( f aiH^{Ai)) n{E + L-\-l, 1)) ^ 1 X ^.,/.^2^(iogL)^ 
^[\aiHUAi))n{E' + L^\-l,l))^^S ) ' 

This decorrelation estimate means that, up to a sub-polynomial error, the 
probability of obtaining simultaneously two eigenvalues near two distinct en- 
ergies is bounded by the product of the estimates given by Wegner estimate 
for each of these two energies. Roughly speaking, two eigenvalues of our 
model near two distinct energies behave like two independent random vari- 
ables. 

Thanks to Theorem 1.2 and an argument in [9], the asymptotic independence 
of the weak limits of S(^, E, u, A) and S(^, E', u, A) associated to the model 
(1.1) is obtained: 

Theorem 1.3. Pick two positive, distinct energies E and E' in the localized 
regime such that iy{E) > and v{E') > 0. 

When |A| — > +oo, the point processes S(^, i?, A), and E' ,u, A) con- 
verge weakly respectively to two independent Poisson processes on M with 
intensity the Lebesgue measure. 

That is, for {Uj)i^j<^j, Uj C M bounded measurable and Uj' Uj = if 
j ^ j' and {kj)i^j<^j G N'^ and (f^j)i^j<j', U'j C M bounded measurable and 
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f/j, n t/j = z/ J ^ j' and (A;;.)i<KJ' ^ N-^ 




/ 


'#{j;0(i?,w,A)Gf/i} =A;r 




P 


< 


#{j;0(i?,a;,A)GC/,} =fcj 
#{j;0(i?',w,A)Gf/(} =A;1 


> 




V 


.#{j;e,(i?',a;,A)Gt/S} =fcS. 


I 



n 



\u. 



k'\ 



converges to as |A| goes to infinity. 

To prove Theorem 1.2 for the model (1.1), we follow the strategy introduced 
in [9]. The key point of the proof of decorrelation estimates for the ID 
discrete Anderson model in [9] is to derive that the gradient with respect to 
u of two eigenvalues E{u) and E'{u) near two distinct energies E and E' 
are not co-linear with a good probability. In the discrete Anderson case, this 
statement will hold true if the gradients of E{u) and E'{uj) are distinct. In 
deed, the gradients of E{uj) and E'{u) have non-negative components and 
their norm are always equal to 1. So, if they are co-linear, they should 
be the same. 

Unfortunately, the norm of the gradient of an eigenvalue of our finite- 
volume operator Hi^{A) is not a constant w.r.t. u anymore (it is even not 
bounded from below by a positive constant uniformly w.r.t. u). Moreover, 
to prove the above key point for the discrete Anderson model, [9] exploits 
the diagonal structure of the potential which can not be used for the present 
case. So, a different approach in the proof is needed to carry out Theorem 
1.2. This approach is contained in Lemma 3.3. 

In addition, according to [8], a result similar to Theorem 1.3 can be obtained 
for our model: 



Theorem 1.4 (Theorem 1.12, [S]). Pick < Eq ^ I such that the density of 
states V is continuous and positive at Eq. 

Consider two sequences of positive energies, say (-E'a)a, {E'^)a such that 



1. Es 



2. \K\\N{E^)-N{E'^ 



> Eq and E', > Eq, 



A^Zd 



-> +00. 



Then, the point processes S(^, E'a, cu. A) and S(^, E'j^, u, A) converges weakly 
respectively to two independent Poisson point processes in M with intensity 
the Lebesgue measure. 

In Theorem 1.4, instead of fixing two distinct energies E and E', one considers 
two sequences of energies {-Ea}, {-E'a} which tend to each other as |A| — >■ oo. 
In addition, we assume that, roughly speaking, there are sufficiently many 
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eigenvalues of i?i^(A) between E/^ and E'j^ as A large. Then, we can derive 
again the asymptotic independence of two point processes associated to E\ 
and E'^. 

Besides, it is known that the existence of an integrated density of states de- 
fined as in ( 1.2) implies that the average distance (mean spacing) between 
eigenlevels is of order |A|~^. 

Thus, Theorem 1.4 shows that, in the localized regime, eigenvalues separated 
by a distance that is asymptotically infinite with respect to the mean spac- 
ing between eigenlevels behave like independent random variables. In other 
words, there are no interactions between distinct eigenvalues, except at a 
very short distance. 

Finally, in Section 5, we will see that a similar statement to Theorem 1.3 
holds true as we consider not only two but any fixed number of distinct en- 
ergies. 

Notation: In the present paper, we use Dirac's notations: \x) = 6x, and 
\x){x\ = {6x, ■)Sx is the projection operator on 6x- 



2 Preliminaries 

We recall here a Wegner-type estimate and a Minami-type estimate for the 
model (1.1) which are essentially important for the proofs of Poisson statistics 
as well as those of decorrelation estimates. 

Theorem 2.1. [Wegner estimate, Theorem 2.1, [6]] 

Fidtst{E,a{H^{A))) ^ e)^^^^^^e\A\ 

hi — e 

for all intervals A C Z and < e < E . 

Roughly speaking, the above Wegner estimate means that the probability 
that we can find an eigenvalue of the finite- volume operator H^{A) in a 
interval (suppose that the length of this interval is small) is bounded by the 
length of this interval times |A|. 
Following is the Minami estimate: 

Theorem 2.2. [Minami estimate. Theorem 3.1, [(>]] There exists C > such 
that, for all intervals J = [a, b] C (0, +oo), and A G 1^, we have 

P(#Mif^(A)) n J} ^ 2) ^ C| JplAp. 

It is not hard to realize a similarity between the Minami estimate and the 
decorrelation estimate in Theorem 1.2: the Minami estimate can be inter- 
preted as a decorrelation estimate for close together eigenvalues and can be 
used to get the counterparts of Theorem 1.1 and 1.3 when E and E' are 
closer and closer to each other as |A| — )■ oo (c.f. [8]). 
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We should remark that although the almost sure spectrum S of the operator 
(1.1) contains the point (the bottom of E), the above Wegner and Minami 
estimates do not work at 0. That is why, in the present paper, we will only 
consider energies who are strictly positive in S. 

We would like to remind readers of the precise definition of the localized 
regime or the region of localization. 

Proposition 2.3. (c.f. [ '], ['>]) Let I be the region of "E where the finite 
volume fractional moment criteria of [.:'] for H^{A) are verified for A suffi- 
ciently large. 

Then, there exists u > such that, for any p > 0, there exists q > and 
Lq > such that, for L ^ Lq, with probability larger than 1 — L"^, if 

1- '^n,u) is a normalized eigenvector of H^^{Ai) associated to an energy 

2- Xn,u) E Al is a maximum of x ^ \^n,ui{,x)\ in Al, 
then, for x G A^, one has 

The point Xn,uj is called a localization center for ipn^uj or En^uj- 

Remark that an eigenvalue of Hi^{A) can have more than one localization 
center. Nevertheless, it is not hard to check that, under assumptions of 
Proposition 2.3, all localization centers associated to one eigenvalue are con- 
tained in a disk of radius ClogL with some constant C. 
Hence, to define the localization center uniquely for an eigenvalue, we can 
order these centers lexicographically and choose the localization center asso- 
ciated to an eigenvalue to be the largest one. 

3 Proof of the decorrelation estimate 

Pick two distinct, positive energies E, E' in /. Let Jl = E + L~^[—l, 1] and 
J'^ = E' + L-^[-l, I] with L large. 

Lemma 3.1. Suppose that u ^ E{uj) is the only eigenvalue o/if^(A) in J^. 
Then 

1. E{uj) is simple and uj ^ E{uj) is real analytic, and if u ^ is the 
corresponding real-valued, normalized eigenvector, it is real analytic in 

CO. 

2. \\V.E{uj)h. = 2 E,eA lin^^f ^here n, = - (7 + 1))(7 " (7 + 1)1 
is a projection in /^(A). 



3 PROOF OF THE DECORRELATION ESTIMATE 



7 



3. Hess^E{uj) = {h^^js)'^^^ where 

• h^^p := -ARe{{H^{A) - E{u)y'^^,^^). 

• ip-y '■= ~ n^y9 = — n(^)± (n^</?) where '^{^p)^- is the or- 
thogonal projection on {(p)^. 

Proof. (1) is true from the standard perturbation theory (see [10]). 
Now we will prove (2). Starting from the eigenequation: 

H^iA)ip = Eioj)ip (3.1) 

we have, for all 7 G A, 

d^^E{uj) = {d^^{H^{A)^),^) + {H^{A)^,d^^^) 

= {d^^iH^{A))ip,ip) + {H^iA)d^^ip,ip) + {H^{A)ip,d^^ip) 
= {d^^{H^{A))^.^) + {d^^^.H^m^) + {H^{A)^^d^.^^) 
= {d^^{H^{A))^,^) + E{uj){{d^^^,^) + {if^d^^^)) 

where the last two equalities come from the symmetry of H^^[A) and (3.1). 
Noting that 

{du;^<P,<p) + {^,d^_^^) = 29^J|v9f = 0. 

Hence, 

d^^E{u) = {d^^iH^{A))ip, if) = 2||n^(^f ^ 0. (3.2) 

Besides, it is easy to check that Il-y = 11* = 11^. Hence, H^ is an orthogonal 
projection. 

From (3.2), we have the following important equality: 

J2^id^,E{uj) = E{u) (3.3) 

which characterize the form of our operator. 

Finally, we give a proof for (3). By differentiating both sides of (3.2) with 
respect to cu-y, we have 

dl^Eiu) = 2(9^^ (H,^),^) + 2(H,^,9^^^) (3.4) 
= 2(H^(9^^(y9, if) + 2(H^(y9, d^^(p) = 4 Re(H^v5, d^^ip) 

Next, we will compute dt^^(p. 

Differentiating both sides of (3.1) with respect to to get 

{d^^H^{A))ip + H^{A)d^^ip = d^^E{uj)ip + E{io)d^^ip 

= 2(H^v^, Lp)(f + E{uj)d^^Lp 

Therefore, 

[H^{A) - E{u)]d^^^ = 2{U^^, - (d^^H^iA))^ = 2(^{U^^, - U^^^ . 
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Observe that := (n^(/?, ip)(p — Il^(p G (v^)"*", and [Hi^{A) — E{uj)] is invertible 
in the subspace (ip)^ of /^(A), we get 

d^^ip = 2{H^{A) - E{u)Y\{Tl^^,p>)^ - n^</.). (3.5) 

From (3.4) and (3.5), we obtain 

dlE{uj) = 4Re(n^<^, {H^{A) - E{uj)Y\{U,^, ^)p> - U,p>)). 

Thanks to (3.5), we have 2(^H^{A) — E{u)) ^ (^{Il^ip,ip)ip — Il^ip) is orthogonal 
to if. We therefore infer that 

dl^Eico) = -ARe{ {H^{A) - E{uj)) " V^7>- 
Repeating this argument, it is not hard to prove that 

dl^^^E{uj) = -4Re((i7.(A) - ^(a;))"'^,, ^/3> 
for all 7, /3. So, we have Lemma 3.1 proved. □ 

Next, we would like to state here a "lower bound" for normalized eigenfunc- 
tions of i^aj(A) in the sense that there exists a large subset J in A such that 
the components {u{k))k£j of a normalized eigenfunction u of Hi^{A) can not 
be too small. 

Lemma 3.2. Let A = [—L,L] be a large cube in Z and (3 G (1/2, 1). Then, 
there exists a point ko in A and a positive constant c such that v?{k) +v?{k + 
1) ^ e"'^'^ for all \k — k^l ^ cL^ when L is large enough. 

Proof. We rewrite the eigenequation corresponding to the eigenvector u and 
eigenvalue E at the point n by means of the transfer matrix 

( u{n) \ 
\u{n - I)) ■ 

Now, we denote the transfer matrix 

UJn + — E —UJ n-l \ 

1 ; 

by T(?T,, E) and denote by v{n) the column vector = {u{n + 1), u{n)y. 
Hence, for all n greater than m we have v{n) = T{n,E) ■ ■ - Tin — m + 
1, E)v{m). Remark that the transfer matrices are invertible. Moreover, they 
and their inverse matrices are uniformly bounded by a constant C depending 
only on the support of ujj, because Uj G [ao, /3o] with < < /^o- Thus, 

\v{n)\ ^ C\"\\v{m) \ = e'°s^l"-'"l|w(m)| = e^l"-'"l|w(m)| 



u{n + r 

u(n) 



1 
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for all n, m in A. 

Assume that |w(/i;o)| is the maximum of Hence, 

|^^(A;o)| ^ ^ 



/2L 

from the fact that YljeA ~ 2. Thus, for any k > 0, we obtain 

V2L 



when \k — ko\ ^ kL^ and L sufficiently large. 
Hence, we have 

u\k) + u\k + 1) ^ e-"^^' 

for all Ifc — A;o| ^ kL^. So, by choosing k = — , we have Lemma 3.2 proved. □ 

V 

The following lemma is the main ingredient of the proof of the decorrelation 
estimate as well as the heart of the present paper: 

Lemma 3.3. Let E, E' he two positive energies in the localized regime and 
/3 > 1/2. Let ci,C2 > and denote by P* the probability of the following 
event (called (*)): 

There exists two simple eigenvalues of Hi^{Al), say E{uj), E'{u) such that 
\E{u) -E\ + \E'{cu) - E'\ ^ e-^^ and 

\\VUclE{u)-clE'{uj))\\,^cle-'^'. 
Then, there exists c > such that 

Remark: There is a slightly but really important difference between the 
above lemma and Lemma 2.4 in [9] where ci = C2 = 1. 

In fact, in the proof of the decorrelation estimate for our model, we will use 
the above lemma with c^,c| are respectively which are two distinct, 

positive numbers. This difference results from the lack of the normalization 
of II V-E'(a;)||i for our model. 

Moreover, we will see in the remark at the end of this section that, for the 
model (1.1), if Ci = C2, P* is exactly equal to when L large enough. 
We will skip for a moment the proof of Lemma 3.3 and recall how to use 
the above lemma to complete the proof of Theorem 1.2. This part can be 
found in [9]. We repeat it here (with tiny modifications) just for the readers 
convenience. 
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Proof. Recall that Jl = E + L-'^[-l, 1] and J'^ = E' + ^-^[-1, 1] with L 
large. 

Let cL°' ^ / ^ L"/c where c is fixed. By the Minami estimate, we have 
P(#{a(i/^(AO) n Ja ^ 2 or #{a(i7^(A0) n J'^} ^ 2) ^ (1/1)'. 

Hence, it is sufficient to prove that Po ^ C(//L)2e(i°g^)'' where 

Po := P(#{a(i/^(A,)) n Jl} = 1; #{a(i7^(A0) n J[} = l). (3.6) 

The crucial idea of proving decorrelation estimates in [9] is to reduce the 
proof of (3.6) to the proof of a similar estimate where A; is replaced by a 
much smaller cube, a cube of side length of order logL. Precisely, one has 
(c.f. Lemma 2.1 and Lemma 2.2 in [9]): 

Po < Cil/Lf + C(//7)Pi 

where / = C log L and 

Pi := F{i^[aiH^{Aj)) n Jl] > 1 and #[a(i/^(A,-)) n J'l] > l) 

where Jl = E + L-i(-2,2) and J' l = E' + L-^{-2,2). 
To complete the prove of Theorem 1.2, one proves that 

Pi ^ C(l/LyJ'. (3.7) 

Combining the Minami estimate and the following lemma 
Lemma 3.4. [Lemma 2.3, [9]] There exists C > such that 

C 

\\Hess^{E{uj))\\i^^ii ^ 

dist[E[(xi), a[H^[A)) \ E{u)) 

we obtain that 

V 1^ \\Hess^{E{u))\\io.^ii ^ e ' } J 
Hence, for e G (4L~^, 1), one has 

Pi ^ CePL-^ + P, (3.8) 

where P, = P(fio(e)) and 

aiH^A^)) n Jl = {Eiu)} 
I {E{u)} = a{H^{A^)) n{E-Ce,E + Ce) 
a{H^{Aj))nJ'L = {E'{uj)} 
{E'{uj)} = a{H^{Ai)) n {E' - Ce, E' + Ce) 
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On the one hand, keeping the same notations as in [9], one puts A 



1 



For 7 and 7' in A, one defines 



nl'J (e) = fio(e) n {CO I \ J,,y{E{u),E\uj))\ ^ A} 

where J.y^.y'{E{uj), E'{uj)) is the Jacobian of the mapping 

{u^,ujy) I-)- {E{uj), E'{uj)). 

Hence, we can dominate Pe by 

P,^ J]P(02;/(e))+P. 
75^7' 

where P^ is the probabihty of the following event 

V:= {uj\ |J^,y(E(u;),E'(a;))| ^ A for all 7,7 G A;}. 
Moreover, from Lemma 2.6 in [9], we know that 

Hence, 

P, ^ Cl^L-^\-^ + ¥r. 
Combining (3.8), (3.9) and choose e := L"^A~^, we infer that 

Pi ^ C(l/LfJ^ + Fr. 



(3.9) 



(3.10) 



Finally, we will use Lemma 3.3 to estimate P,.. 

To do so, for each u & V, we rewrite the Jacobian J^^y{E{u), E'{u)) as 
follows: 



J^,y{E{uj),E'{uj)) 



d^.E{u) d^_^,E{u) 
d,.,E'(u) a, ,E'(u). 



E{u)E'{u) 



if E{u) and E'{uj) are non-zero. 

Observe that the sum over 7 or 7' of each entry in the second line of the 
above 2x2 matrix is equal to 1 by (3.3), we can use the following lemma: 

Lemma 3.5 (Lemma 2.5, [9]). Picku,v G (M"*")^" such that \\u\\i = \\v\\i = 1. 
Then, 



max 



Uj Uk 
Vj Vk 



1 



\u — v\ 
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to deduce 

\\VU^Eiu)-^E\co))\U^e-^' 

provided that 

\E{uj) -E\ + \E'{uj) - E'\ ^ e-'' 
and / = O(logL) is sufficiently large. 

Thus, thanks to Lemma 3.3, for L sufficiently large, and with the remark 
that / ^ clogL, (3 > 1/2, we have 

which completes the proof of Theorem 1.2. □ 

To complete the present section, we state here the proof of Lemma 3.3: 

Proof of Lemma 3.3. Let u H- u{uj) and u H- v{uj) be normalized eigenvec- 
tors associated to E{u!) and E'{uj). By Lemma 3.1, we have 

W^E{uj) = (2||n,«f )^^^ and V^E'{u) = (2||n,^f )^^^ 

Now, We introduce the linear operator T from /^(A) to /^(A) defined as 
follows 

Tu{n) = u{n) — u{n + 1) 

where u = {u{n))n G /^(A). Recall that A = A^ = Z/LZ, i.e. we use periodic 
boundary conditions here. 
Assume that 

c?e-^' > WVUclEiu) - clE'iu))\\i = \iciTuin)f - {c2Tv{n)f\ 

n 

= \ciTu{n) — C2Tv{n)\\ciTu{n) + C2Tv{n)\. 

n 

Hence, e'^^ ^ \Tu{n) - cTv{n)\\Tu{n) + cTv{n)\ with c = cs/ci > 0. 
Then, there exists a partition of A = {—L, . . . , L}, say P C A and Q C A 
such that P U Q = A, P n Q = and 

• for n eV, \Tu{n) — cTv{n)\ ^ e~^^ ^'^ 

• forn e Q, \Tu{n) + cTv{n)\ ^ e~^^^^ 

In the rest of this proof, we put v{n) := cv{n). This abuse of notation changes 
nothing thanks to the linearity of the operator T. 

From Lemma 3.2, there exists a Ci > and an interval J of the length ciL^ 
which satisfies 

|m(A;)|2+|u(A; + 1)P ^e-^'/^ (3.11) 
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for all k G J. 

Now, we want to study the properties of the above partition in the interval 
J. 

We assume that 

p n J = UPj and Q n J = UQj 
where Vj and Qj are intervals in Z. 

Step 1: We'll show that the length of any Vj or Qj is inferior or equal to 3. 
Indeed, assume by the contradiction that there exists an interval Vj = {n — 
2,n — l,n,n + l,n + 2 . . . m} with \Vj\ = \m — n + 2\ ^4. So, m ^ n + 2. 
First of all, the following relations are true for the points in Vj 

Tu{n - 2) = Tv{n - 2) + 0{e~^^/^). 
Tu{n - 1) = Tv{n - 1) + 0{e~^^''^). 

Tu{n) = Tv{n) + 0{e~^^/^). (3.12) 
Tu{n + 1) = Tv{n + 1) + 0{e-^^/^). 
Tu{n + 2) = Tv{n + 2) + 0(6-^'"/^). 

Combining the eigenequations for u and v at the point n — 1 

Eu{n - 1) = ujn-iTu{n - 1) - uJn^2Tu{n - 2) + 0(6"^^/^) (3.13) 

E'v{n - 1) = a;„_iTt;(n - 1) - cUn-2Tv{n - 2) + 0(6-^"/^) 
and relations in 3.12, we obtain 

Eu{n - 1) = E'v{n - 1) + 0(6-^"/^). (3.14) 

Similarly, we have 

Euin) = E'vin) + ©(e'^'/^). (3.15) 
Combining (3.14), (3.15) and the second relation in (3.12), we obtain 







(i -!)"(")= 









which implies that 

Tu(n- l)^Ce-^'/2 (3.16) 
where C is a positive constant depending only on E, E' and the support of 

UJ. 

In other words, we showed that, if we have three consecutive points in some 
interval Vj (n — 2, n — 1, n in the above case), then the middle point of those 
three points always satisfies the inequality (3.16). 
Repeating the above argument for other points in Vj, we obtain that 

TM(r^)^Ce-^'/^ (3.17) 
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and 

TM(n+l)^Ce"^'/^ (3.18) 

Therefore, combining the equations (3.16), (3.17), (3.18) and the eigenequa- 
tions for u at the point n and n + 1 

Eu{n) = uJnTu{n) — uJn-iTu{n — 1) + 0{e~^'^^'^), 
Eu{n + 1) = uJn+iTu{n + 1) - uJnTu{n) + 0{e~^^''^), 

we infer that there exists a positive constant C being independent of L such 
that 

|M(n)p + |M(n + l)p^Ce-^'' 
which contradicts (3.11) if we choose L large enough. 

Hence, an interval Vj or Qj can not contain more than four points in Z and 
we have the statement in Step 1 proved. 

Step 2: Consider four arbitrarily consecutive points in J, say n — 2,n — 
1, n, n + 1. The following are all possible configurations for these points with 
the remark of the equivalent role of Vj and Qj. 

First case: We consider the case when three of these points belong to an 
interval Vj and the other belongs to Qj , say n — 2,n — l,n G Vj and 
n + lE Qj. 

First of all, we have the same eigenequations as in (3.12) with respect to 
n — 2,n — l,n and the following 

Tu{n + 1) = -Tv{n + 1) + 0(6-^"/^). 

Secondly, we have two eigenequations for u and v at the point n + 1 

Eu{n + 1) = Un+iTu{n + 1) - ujnTu{n) + 0(6"^"/^) (3.19) 

and 

E'v{n + 1) = Un+iTv{n + 1) - UnTv{n) + ©(e"^"/') 
which can be rewritten as follows 

E'v{n + 1) = -cOn+iTu{n + 1) - uJnTu{n) + 0{e-^^/^). (3.20) 

Thirdly, repeating the argument of Step 1, we can obtain (3.14) and (3.15) 
again. So, combining with two eigenequations for u at the point n — 1 and 
n, we will have the following system of linear equations (10 equations, 10 
variables) : 

AoU = bo 

where U = {u{n — 2), u{n — 1), ... , u{n + 2), v{n — 2), v{n — 1), ... , v{n + 2)), 
^0 = (^o)i=sj<io with ||6o|| ^ ce~^^ and Aq is the following 10 x 10 matrix: 
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OJr. 
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-1 
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-1 








o\ 
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-1 
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-1 
















-1 


1 



























-1 




































E' 



























-1 








0/ 



Let Aq be the comatrix of Aq and po be max{l, HAqH}. 
We will prove that | detAo] ^ poe~^ 
Indeed, assume by the contradiction that | detAo] > poe' 
Therefore, Aq has the inverse satisfying 

Po 



.L'3/4 



1^0 I 



det An 



^ e 



Hence, 



1^0 I 



Poe 



But this is impossible according to (3.11). Consequently, | det Aq\ 
Remark that po only depends on the support of Uj and E,E'. 
Now, we would like to compute explicitly the determinant of the matrix Aq 
which is a polynomial of at most four variables uj with j = n — 2,n + 1. 
Then, combining with the fact that this determinant is small, we infer some 
restrictive conditions on u. 

Put Aq = {ttij), we'll give here some details of computing the determinant of 
the matrix Aq by hand: 

(A mathematical software like Maple or Mathematica might be useful for 
checking the final result of this computation). 

First, expand this determinant by its sixth and last column and then by its 
first column to get 

I det Aq I = Un-2 1 det Bq \ 
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where Bq is the 7x7 matrix defined by 



1 


-1 








1 


-1 


o\ 





1 


-1 








1 


-1 




l^n-l + U^n — E 



















E 











-1 










iOn + UJn+l — E 




























E' 


E 
~E' 











-1 





V 



V 

Now, we compute the determinant of Bq. 

Take the sixth row minus the fifth row and take the first row plus the last 
row. Next, multiply the second row by E' and take it plus the sixth row. 
Finally, expand the determinant of Bq by the forth, the fifth and the last 
column to get 

I det Ao\ = —ujn^2UJn+i\ det Co\ 
where Co is the 4x4 matrix defined by 



/ E 




E 



-1 

E' + 2ujn E-E' - 2uJn 

E -Un 





-1 

E' 


-1 



Finally, by an explicit computation for the determinant of Co, we obtain that 



AE 

det Aol = —{E + E')uJn-2UJn+l 



E' -E 



4 



Poe 



Therefore, with the remark that E, E' are positive and {u)j}jizz are bounded 
from below by a positive constant, the following condition on uj holds true 
when L sufficiently large: 



E' -E 



(I) 



Second case: Now assume that n — 2, ?i — 1 G Qj and n,n + 1 G P-,-. 
Repeating the argument in First case, we have a 10 x 10 system of linear 
equations 

AiU = bi 

where U = {u{n — 2),u{n — 1), ... , u{n + 2), t>(n — 2), w(n — 1), ... , v{n + 2)), 
bi = (^)i^j^io with ^ lOe"^''/^ and Ai is the following 10 x 10 matrix: 
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-1 







— LOyi—I 


OJn-l + OJn ~ E 


— 














UJyi — ^n— 1 


— 















—UJn 


OJn \ CUn-l~1 


E - 


'-^71+ 1 













E 

'W' 

























n 


E 

~E' 


n 


n 




n 










-1 1 








\ 








-1 


1 
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-1 
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-1 





























-E' 






































-1 





























-1 








0/ 



Denote by Ai the comatrix of Ai and put pi := max{l, As in First 

case, we obtain that 

I detail ^ pie-^^/\ 
The determinant of Ai can be computed as follows: 

First, expand this determinant by its sixth and last column and then by its 
fifth and first column to get 

I det Ai I = Un-2^n+i I det Bi \ 

where Bi is the 6x6 matrix defined by 



/ 1 


-1 





-1 


1 


o\ 





1 


-1 





-1 


1 




Wn-l + OJn — E 


-OJn 













OJn — OJn-1 


-OJn 





-E' 











E 
E' 








-1 


E 

\ ~E^ 








-1 





V 



Second, take the first row of matrix B minus its last row and take the second 
row minus the fifth row, then expand the determinant of Bi by its forth and 
sixth column to obtain 



detail 



OJr, 



-20Jn+i \ detCil 
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where Ci is the following 4x4 matrix 

-1 



^-1 



1 





E 

-1 

E' 



1 

-E' 
-E' 



Finally, by an explicit computation for the determinant of the matrix Ci, we 
infer that 



detail 



4E 



{E + E 



•l\2 



4 



^ Ce 



-LP /A 



Hence, take into account that Uj ^ ao > for all j & 1^ and E, E' > 0, we 
have the following condition on the random variables: 



{E + E 



A2 



4 



^ Ce 



{II) 



as L sufficiently large. 

Third case: Suppose that n — 2, n — 1 G Qj, n G Vj and n + 1 G Qj+i- 
Proceeding as in the previous cases leads us to consider the following system 
of linear equations 

A2U = 62 

where A2 is the 10 x 10 matrix defined by 



/ 1 


-1 
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-1 
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-1 
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-1 





—^n-l 


+ UJn — E 

























—Un 


LOn + tOn+1 " E 












(jJn — UJn+1 






W„,-2 + — E 

E 
~E' 










V ° 












-1 

























-1 
1 




-E' 


















o\ 




1 







-E' 

0/ 
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and \detA2\^Ce-^^/\ 

Then, by an explicit computation, we obtain that 

E' -E 



I det A2I = 4:E{E + E')ujn-2 
which immediately yields that 



4 



^ Ce 



-LP/4 



E' -E 



4 



^ Ce 



{III) 



as L — > +00. 



Forth case: Suppose that n — 2 G Qj, n - 


- 1 e -Pj, 


n G 




dliU 


n + 1 G Vj+i. Then, U = {u{n — 2), u{n — 1), . 


. . , u{n + 2 


),v{n 


-2), 




1), . . . , f (n + 2)) satisfies the following system of linear equations: 




AiU = b-^ 










where ^3 is the 10 x 10 matrix defined by 










/I -1 












1 -1 
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-1 




— t^n-l + UJn — E 











—OJn-l OJn-l — 


UJn 









—Un 00 


n + — 


E - 






UJn 










—Un~2 ^n~2 + ^n~l ~ E —UJn-l 


























-1 1 
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-1 













-1 
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-1 





















-E' 






























-E' 





















-E' 








0/ 



and I det A3 1 ^Ce-^''/^ 
We compute 

I det A3 1 = EE'Un-2 X Un+l X 



icon-i + E' - E 



E' -E + iujn 



Hence, at least one of the two following conditions on u must be satisfied: 
E' -E 



4 



^ Ce-^'/^ {IV) 
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E' -E 

Un + 



4 

Fifth case: Assume that n — 2, n — 1, n, + 1 G Vj. Hence, by Step 1, 
n + 2 G Qj. We consider four points n — l,n, n + l,n + 2 and come back to 
First case. 

Last case: Suppose that n — 2 G Qj, n — l,n & Vj, and n + 1 G Qj+i- 
We consider the point n + 2. 

If 77, + 2 belongs to Qj+i, we consider four points n — l,n,n + l,n + 2 and 
come back to Second case (by the equivalent role of V and Q). Otherwise 
n+2 belongs to Vj+i. In this case, we consider four points ra — 1, n, n + 1, n + 2 
and come back to Third case. 

To conclude, the random variables Uj must satisfy at least cL'^ conditions 
of one of the types (I)-(IV). From the fact that Un are i.i.d. and possess a 
bounded density, the conditions (I)-(IV) imply that the event (*) can occur 
for a given partition V and Q with a probability at most, (e~'^^'^)^'^ = e~^^^^ , 
for some c > 0. Hence, 

as the number of partitions is bounded by 2^ and (3 > 1/2. 

We therefore have Lemma 3.3 proved. □ 

Remark: Thanks to the equality (3.3), it is not hard to derive the following 
estimate for the model (1.1), 

A E 

-^|Ar'/' ^ W'^UEico) - E\u))h ^ WVUEiu) - E'iu))h (3.21) 

provided that \E{u]) — E\ + \E'(lo) — E'\ ^ e~^^ and AE is the distance 
between E and E'. 

The above estimate reads that the /^—distance of the gradients of E{oj) and 
E'{(jj) is bounded from below by a positive term that is not exponentially 
small w.r.t. the length of the interval A. 
Now, let Ci = C2 in Lemma 3.3. 

Under the hypotheses in Lemma 3.3, the estimate (3.21) and Lemma 3.5 
immediately imply that 

that is impossible if |A| is sufficiently large. Hence, for ci = C2, P* is equal 
to . 

Finally, we would like to note that an estimate as (3.21) for the discrete 
Anderson model only holds true for two distinct energies sufficiently far apart 
from each other. Moreover, that kind of estimate enable us to prove the 
decorrelation estimate for the discrete Anderson model in any dimension 
(c.f. Lemma 2.4 in [9]). But it is not the case for the model (1.1). 
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4 Comment on the lower bound of the r.v.'s 

In this section, we want to discuss how to relax the hypothesis on the lower 
bound of random variables {uj}j^z- 

Assume that all r.v.'s Uj are only non-negative instead of being bounded 
from below by a positive constant. 

Thanks to the Taylor expansion, the common distribution function F{t) of 
{ujj}j^z always satisfies the following equality: 

= p(0)t + o(t) (4.1) 

near where p is the density function of the r.v.'s {wjjjgz- 
Hence, 

m > (4.2) 

near provided that p(0) is positive. 

We will show that decorrelation estimates can be obtained by using (4.2) i.e., 
by using the assumption p(0) > 0. 

Pick A = [—L, L] as usual and f3 a fixed number in (0, 1) which will be chosen 
later. Then, we have 2L + 1 — sub-intervals of length in A. 
Using the independence of the r.v.'s Uj, we infer that 

:= P(3 an interval J of length in A such that ^ e"^'"^^)' "^J ^ J ) 
^ (2L + 1 - L^)(l - P(a;o < e~^^"^^^'))^'' ^ 2L(1 - F{e~^^°^^^')f' 
where 6 is an arbitrary, fixed number in (0, 1/3). 

Next, we put x := —F{e~^^°^^^ ) G (—1,0), and rewrite the quantity (1 — 
^(g-(iogL)^))L'9 in the form {{1 + xy/'^Y^^ 

We observe that the function f{x) = (1 + xY^^ is decrease in (—1,0) and 
limrc^of^x) = e, hence f{x) ^ e Vx G (—1,0). 
Combining this with (4.2), for L large enough, we get 



Pi ^ 2Le-^'^(^"'°^'^') ^ 2Le-^»^''^""°^''' < 2Le-^^'^ < e-'"^'' (4.3) 

where Cq,c,/3i are positive constants and Pi is smaller than /3. 

In other words, with a probability 1 — e~'^^ ^ , there exists a sub-interval J of 

side length L'^ of A such that 

^ e~(^°si)^ vj G J. 

We can therefore use the same argument in the proof of Lemma 3.2 to get: 

Lemma 4.1. Let A = [—L,L] be a large cube in Z and e be a fixed number 
in the interval (0,/3). 

Then, with the probability l—e~'^^'^^ , there exists a point ko in A and a positive 
constant such that v?{k) + v?{k + 1) ^ e~^^ for all\k — k^l ^ c^L^^^''"'' as 
L is large enough. 
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Hence, the upper bound for the probabihty P in Lemma 3.3 is now B 
c2^e~'^^^ '^'^^^ if we choose e to be /3/2 in Lemma 4.1. 



V 17 - 1 

Observe that /3^ + /3/2 > 1 provided that 1 > (3 > . With this choice 

of P, the upper bound B is exponentially small as L tends to plus infinity. 
Then we can proceed again as in Lemma 3.3 and remark that the error term 
g-cL/^i ^ (//L)^e^'°^^^'', we obtain again a decorrelation estimate for the 
model (1.1): 

Theorem 4.1. For (3 G {-~^ — ij, a G (0, 1) and I = 0{L°'), we have 

^l\a(i/.(A,))n(E' + L-i(-l,l))^0/j ^^^^/^^ " -o(//L). 
when L is sufficiently large. 



5 More than two distinct energies 

In this section, we would like to show that, following an argument in [9], 
we can use Theorem 1.2 to prove the asymptotic independence for any fixed 
number of point processes, i.e., for a fixed number n ^ 2 and a finite se- 
quence of fixed energies {-Ei}i<i<n in the localized regime, n point processes 
S(^, Ei, uj, A) defined as in (1.3) converge to n independent Poisson processes. 
For example, we will prove in detail the case of ra = 3 with three distinct 
energies E,E', E" . 

Consider non-empty compact intervals (f/j)i^j<j, {U'j)i^j^j',{Uj)i^j<^j" and 
integers (/^jOisr^j, {kj)iiiji:j', (^j)i=£jXJ" as in Theorem 1.3. 
Using notations in [9], we pick L and / such that (2L + 1) = (2/ + 1)(2/' + 1), 
cL" ^ Z ^ L"/c for some a e (0, 1) and c> 0. 

We decompose the interval A := [—L,L] = IJ|7|^z'^K7) where A;(7) := 
(2/ + l)7 + Ai. 

Next, for A' C A and f/ C M, consider the r.v.'s 

{1 if H^{A') has at least one eigenvalue 
in E + {u{E)\ A' \)-'U, (5.1) 
otherwise 

and put J:{E,U,1) := j:^^^^, X{E,U, Ai{-f)). 

First of all, [Lemma 3.2, [9]] is the first ingredient of the proof which tell us 
that we can actually reduce our problem to consider eigenvalues associated 
to smaller intervals. This lemma is still true in the three -energy-case. 
Hence, to complete the proof of the stochastic independence w.r.t. three 
processes, the only thing left to do is show that the joint finite dimensional 
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distribution function of three point processes tends to the product of three 
finite dimensional distributions associated to each of three energies, i.e., 




S(E,f/i,/) = fci, 

u- j:{E',u[,i) = k[, 

T.{E\Ull) = kl 



,j:{E,u{j),l) = kj 
,E{E',U'j,J) = k'j, 

j:{e",u'J„,i) = k'}„ 



should be well approximated by the product 



P 



v 



S(^,t/i,/) = 


ki' 












1) 


X P 


(1 


E{E,uU)J) = 


kj^ 









j:{e',u[,i): 

E{E',U'j„l) 



X P 



^ E{E",U'j.,l)= k"j„] ) 



(5.2) 



as L goes to infinity. 

By a standard criterion of the convergence of point processes (c.f. e.g. The- 
orem 11.1. VIII, [3]), the above statement holds true if the following quantity 
vanishes for all real numbers tj,t'j,,t'j„ : 



as L goes to infinity. 

Remark that, from the fact that {A(7)}|^|5;;/ are pairwise disjoint intervals, 
operators {H^{K{^))}\^\<^ii are independent operator-valued r.v.'s. 
We thus have: 



Our goal is to approximate terms of the following form 
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by the product 

J J' J" 



tiX{E',m,a) TT j^^-t';„x{E",u'^„a) 



i j' j" 

as L large enough with the remark that these r.v.'s are not independent. 
To illustrate our computation for that, we will consider just three arbitrary 
Bernoulli r.v.'s Xi, X2, and we will compute explicitly A := E(e^»=i "*"^'). 

i<j ^ ^ i=l 

First of all, we rewrite 

p(.Y. = x. = x3 = o) = i-$:( _^;!o=J'^,) 

and obtain that 

+ J2(<''-*'- - i)p ( jflTofi J ^ j) + {'^'•*"*" - mf]x,). 

i<j ^ ^ i=l 

Next, use 



i=l 



to get 



^ = 1 + J](e- - 1)P(X, = 1) + J](e- - l){e^^ - 1)P [^^^ ^ 
+ I e^Li'^. _ 1 _ ^(e»« _ 1) j P(p|Xi). 



i=l / i=l 
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X = X = 1 

UsinET a similar expansion for all terms of the form P ( * „ / , . ' . 



we 



obtain the following expression of E ^e^^=i"*"^*j 



A = l + J](e°' - l)F{Xi = 1) + J^^e""^ ~ ^)(^''' ~ = = 

3 / ^ 

+ J](e'^»-l)P f|X, 



i=l v=i 
On the other hand, from the observation that 



Ee%Aj = 1 + (e^. _ l)]p{Xj = 1), Vj = {1, 2, 3}, 
one multiplies the three equalities above to get 

3 

■QE(ga,x,^) _ ^ ^ J^^e"' - l)P(Xi = 1) 

i=l 

3 3 

+ ^(e'^' - l)(e'^^ - 1)P(X, = 1)P(X, = 1) + Hie'' - 1) H^l^^ = 1^ 



i<j 



i=l 



1=1 



Hence, the difference between the expectation of e^?=i°'^' and the product 
of expectations e"''^^ is: 

3 
i=l 

^(e«« - l)(e'^^- - 1) [F{Xi = X, = 1) - P(X, = 1)P(X, = 1)] 



i<j 



i=l 



1=1 



Thanks to the decorrelation estimate and the Wegner estimate, we obtain 

P(X, = X, = 1) + P(X, = 1)P(X,- = 1) 
^ C(//L)2(e('°s^)' + 1) ^ C{1/LY+p 



and 



P(f|X,) + P(Xi = 1)P(X2 = 1)P(X3 = 1) 
i=l 



5 MORE THAN TWO DISTINCT ENERGIES 



26 



where p G (0, 1). 

To sum up, we get that: 



i=l 



In other words, E(e^^=i '^'"^*) is well approximated by 11^=1 ^(^"*"^0 with a 
small error term of the magnitude It is easy to see that we can 

obtain the same result for an arbitrarily finite number of Bernoulli r.v.'s. 
Hence, for each I7I ^ /', we have: 

E ^g-E, tjX{E,u{j)n)-~j:j> t'^,x(E\m,rt)~T.j» t'^„x(E",u'^„,'r) 



X (l + O(//L)^+0. 

Besides, we have similar formulas for the expectation of e~^^^^^^^''^^-^^''^\ 
^-j:,t,xiE',u;,,) ^-j:,t,xiE",u^,,)_ 

We therefore infer that: 

E (f^-T.jtjX(EMj),7)-T,j>t'^,xiE\u'^,a)-T.j'>t';,,xiE'',u:.',,,-/) 



= Ee- tjXiE,uU),y)^e- E, t,x(E',u;a)^^- E, t,x{E",u'/a) 
X {1 + 0{1/LY+p). 

An important observation is that we only have /' = 0{L/l) sub-intervals of 
the side length / of the interval A. Hence, as multiplying all above equalities 
side by side over 7, we obtain that 

is equal to the product of 

and an error term of the form 

which tends to 1 as L goes to infinity with x = 0{l/L). 

Hence, the stochastic independence for three point processes w.r.t. three 

distinct energies is proved. 
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